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Abstract. The dephasing and relaxation kinetics of bosonic excitons interacting with a thermal bath of 
acoustic phonons is studied after coherent pulse excitation. The kinetics of the induced excitonic polariza- 
tion is calculated within Markovian equations both for subcritical and supercritical excitation with respect 
to a Bose-Einstein condensation (BEC). For excited densities n below the critical density n c , an exponen- 
tial polarization decay is obtained, which is characterized by a dephasing rate F — jr. This dephasing 
rate due to phonon scattering shows a pronounced exciton-density dependence in the vicinity of the phase 
transition. It is well described by the power law r cx (n — n c ) 2 that can be understood by linearization of 
the equations around the equilibrium solution. Above the critical density we get a non-exponential relax- 
ation to the final condensate value p° with \p(t) \ — \p°\ cx i that holds for all densities. 
Furthermore we include the full self-consistent Hartree-Fock-Bogoliubov (HFB) terms due to the exciton- 
exciton interaction and the kinetics of the anomalous functions Fk = {a^a-h)- The collision terms are 
analyzed and an approximation is used which is consistent with the existence of BEC. The inclusion of the 
coherent x-x interaction does not change the dephasing laws. The anomalous function Fk exhibits a clear 
threshold behaviour at the critical density. 
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1 Introduction 

Bose-Einstein condensation (BEC) is a fascinating topic 
attracting new interest, due to recent experimental obser- 
vation for BEC of atomic ensembles |Q,|2j and evidence for 
BEC of excitons in semiconductors ||[I|,|5|,||]. In spite of 
the many similarities between bosonic atoms and excitons 
in semiconductors with large binding energy, the BEC of 
the latter has the interesting specific feature that the or- 
der parameter of the excitonic condensate is identical to 
its optical polarization. Therefore it is directly accessible 
by relatively simple optical measurements. This fact stim- 
ulated theoretical predictions of unusual nonlinear optical 
properties in the condensed regime Q] . 
In this paper the decay of the polarization p (=order pa- 
rameter) is studied after excitation with a coherent laser 
beam, if one approaches the critical density n c for a BEC, 
where an unusual dephasing kinetics can be expected. Un- 
like most studies of the condensation kinetics |^,^,|l^], here 
a large macroscopic occupation of the condensate is cre- 
ated by the excitation process via allowed optical transi- 
tions. This induced exciton amplitude may then decay to 
its zero or nonzero stationary value. The final stationary 
solution is perfectly described by the equilibrium theory 
of the free Bose gas. 

One of the most promising candidates for an excitonic 
BEC is Cu-2.0 with its extremely stable excitons. Due 
to its dipole forbidden exciton transition, an experiment 
corresponding to our simulations would use two-photon 
transitions for the excitation of ortho-excitons and a time- 
resolved polarization measurement via the same mecha- 



of excitons in the vicinity of n c 

nism |l3| . Although for a detailed quantitative simulation 
of such an experiment, additional effects like polariton ef- 
fects |l2|, lifetime effects jlO| and the subtle excitation 
mechanism will play a role even in the low density limit, 
our main predictions should be observable at least quali- 
tatively also in CU2O close to a BEC. 

2 Incoherent relaxation and dephasing 
kinetics for excitons interacting with a bath 
of acoustic phonons 

ls-excitons with center-of-mass momentum k are treated 
in the boson approximation, i.e. the exciton operators ful- 
fill the commutation relation [ak)«k']_ = ^k,k'- This ap- 
proximation is only justified in semiconductors with large 
exciton binding energy, when naj] -C 1 and mainly ls- 
excitons are excited. These excitons interact with longitu- 
dinal acoustic phonons and a coherent classical light pulse 
via a dipole interaction. The Hamiltonian is given by 

+ (1) 

k q 

+ 7v E^ a k+q a k(6 q + blj - W(dE(t)al + h.c.) 

k,q 

with the dispersion of the excitons and phonons 
h 2 k 2 

ek = ~2^T ' Wq = c ' q ' ' ^ 

respectively, where m is the translational exciton mass and 
c the sound velocity. The deformation potential matrix 
element 

,9 q = G^hujq (3) 

in the long- wave length limit. The interaction constant G 
is given in terms of the deformation potential D and the 
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crystal density p by G 2 = j^. The finite photon momen- 3 Laser pulse induced Bose Einstein 
tum has been neglected. condensation 

We define further the cxciton distribution function = 
(ajjafe) with k and the exciton polarization amplitude 

One gets a closed set of kinetic equations for the po- 
larization p and the exciton occupation n^, by extending 
the standard method for deriving a semiclassical Boltz- 
mann equations also to the order parameter p The 
Heisenberg equations for (ao) and (a^ak) are iterated to 
second order in the interaction potential g q . The higher or- 
der mean-values are factorized and the only correlations 
kept are p and n^, e.g. 



The strength of the laser pulse can be varied to excite den- 
sities below (see Fig.[j]) and above (see Fig. ||) the critical 
density for the BEC. The laser pulse is tuned to the low- 
est ls-exciton resonance and has a duration of 2.5ns. As 
an example ZnSe parameters have been chosen at a low 
temperature of 0.5K. Although the BEC in ZnSe may be 
unlikely, this semiconductor with its dipole allowed exci- 
tons serves well to demonstrate our ideas, which may hold 
for different experimental setups and materials. 



{aa+ab+b) » (1 + {a+ a)) (a) {b+ b) + {a){a + a){b+b). 

Neglecting principal value contributions as well as finite 
lifetime effects and performing the Markov limit, one ar- 
rives at B 



In Fig. [j] and Fig. || a light pulse induced condensa- 
tion is obtained. At excitation densities below the critical 
density the polarization \p\ decays exponentially until it 

vanishes as expected from usual optical experiments (Fig. 
d 1 

= —y 2_j {^ / kk' n k(l + «-k') — (k k')} Q). At stronger excitations above the critical density, the 

k' 

— [Wkon-k — Wok(l + n-k)] \p\ 2 (4) exciton polarization approaches non-exponentially its fi- 



df 2V 

k ' 



® 1 [tjt T ,7 /-, , m nite stationary value (Fig. @) , which is non-zero due to 

-foP= ^72^^ k '° nk ' ~~ Wok'(l + «k')JP J lj^ 

the onset of condensation. This stationary polarization is 

given in terms of the excited density n and the temper- 
The transition rates are given by Fermi's golden rule ature T in thc form of | p 0| = ^ n _ n ^ as expected 



7^ dE o(t) • (5) 



Wkk' = ^ |ek - e k ' 

+ (l + JV k _ k /)tf(e k -e k /-n«k_ k /)l , (6) 



Ak'-k<5(ek — ek' + ?ii-L>k'-k) form the equilibrium theory of the Bose gas. Obviously 

the finite lifetime of the excitons due to spontaneous re- 
combination will eventually cause an additional dephasing 
and a decay of the exciton density. 



where Nk = 0h^ k _ 1 is the thermal distribution of the 



phonons. The delta- functions are broadened into Lorentz 
-ians according to |8j. In what follows (except in the last 
chapter) we will study the consequence of these equations. 
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that the contribution of the phonons to the dephasing rate 
is approximately density- independent. We can reproduce 
this result for moderate densities but do get a considerable 
reduction of the dephasing rate near the critical density as 
shown in Fig.|| The density dependence of the dephasing 
rate can be fitted extremely well by a quadratic law 



r oc (n — n c ) z 



(7) 



Fig. 1. Complete exponential dephasing for excitation with 

n < n c 



which vanishes at the critical density, as also shown in 
Fig. [3| This characteristic slowing-down of the exciton- 
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Fig. 2. Non-exponential approach of the exciton amplitude to 
condensation for excitations with n > n c . Full line - calculated 

polarization, dashed line - value predicted from equilibrium Fi S- 3. Density dependence of the dephasing rate T for n c = 
theory of the Bose-gas 2 ■ 10 cm 

4 Critical slowing-down of the dephasing for 



approaching the critical density from below 

Below the critical density the decay of the polarization is 



dephasing kinetics can serve as a first experimental sign 
for the approach to the BEC. 

A similar result has been observed also for exciton-exciton 
exponential, therefore a dephasing time T 2 and a dephas- scattering within similar rate equations |u| . 
ing rate r = ^ can be defined. We study the exciton- The quadratic law Eq.(^) can be derived analytically 
density dependence of these quantities as one approaches by linearizing Eqs. |i|and fj around the stationary solutions 
the critical density from below. It is generally assumed \p°\ and For the deviation 5p from the equilibrium 
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value p° = one gets: 

^P=^EM' " Wbk'(l + n£0] «P. (8) 



of excitons in the vicinity of n, 
98 



2V 



Therefore one can immediately conclude that the final de- 
cay for n < n c is exponential and furthermore 

1 



r = 4 = ^(i-e^)5> k0 - 



(9) 



On the other hand (1 — e' 3 ' 1 ) oc [n — n c ) 2 holds as an exact 
thermodynamic relationship between the leading terms for 
fi « and n sa n c . However, the relation is relatively 
well obeyed even for large departures from the critical 
values, as can be checked by a simple numerical calcu- 
lation Furthermore, due to the structure of the transition 
rates Eq.(^J), the main contribution in the sum of Eq.(^|) 
stems from the fc-values away from the origin (the non- 
zero solution of the delta- function) . Therefore, for small 
temperatures the ^(-dependence of the sum is weak and 
does not affect the (n — n c ) 2 behaviour of its prefactor. 

5 Power law relaxation to the equilibrium 
condensate for n > n c 

While the polarization \p\ vanishes exponentially at sub- 
critical excitation, it reaches its stationary value \p°\ = 



\Jn — n c 0, when the excited density n exceeds n c . But 
in contrast to the subcritical behaviour, this approach to 
equilibrium of the ideal Bose gas is not exponential any 
more. In contrast the dephasing kinetics is well described 
by a simple power law, 
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Fig. 4. Asymptotic kinetics of \p(t)\ for excitations with n > 
n c , power law fit |5p(t)| oc dotted line. Exponential fit 
dashed line. 

as shown in Fig. ||. 

The slow non-exponential approach to equilibrium for 
n > n c would also give in experiments sensitive to coher- 
ence a clear signature for the BEC, which may be observ- 
able even when the final condensate is small and cannot 
be detected clearly. 

Above the critical density the same linearization pro- 
cedure like the one used in the last section but with /j, = 0, 
p = p° + Sp and p° 7^ yields 



(11) 



As shown in Ref. g, the long-time behaviour for dnu is 
not exponential. The spectrum of the linear operator con- 
trolling the t — > oo asymptotics of 8nu is continuous and 
starts from zero. Therefore, it is expected that a power-law 
decay is obeyed, which is in accordance with our numerical 
results. 



\p(t)\=\p \+a/(t + b), 



(10) 
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6 Combination of phonon-scattering with 
HFB-correlations 

The model used in the last sections may be useful to de- 
scribe the incoherent dephasing and relaxation phenom- 
ena of the light-induced BEC. However, it does not ac- 
count for a variety of interesting and important phenom- 
ena, where exciton-exciton correlations play an impor- 
tant role. These excitonic correlations are known to be 
important for the nonlinear optical response of semicon- 
ductors Jl^,|l^,|l^|. Additionally, essential aspects of the 
BEC can be only understood in terms of an interacting 
Bosc-gas. These phenomena are e.g. superfluidity ||,f|], 
the renormalized spectrum of the elementary excitations 
and some interference phenomena based on the Gross- 
Pitajevski equation. These effects are taken into account 
by the self-consistent Hartree-Fock-Bogoliubov treatment 
of the boson-boson interaction that is also called Girardeau- 
Arnowitt approximation. The impact of this approxima- 
tion on the BEC has been reviewed e.g. in Refs. jL8|. We 
will now give a guideline how to include the HFB approx- 
imation in the phonon scattering kinetics in such a way, 
that the results of the last sections are conserved, but 
moreover coherent aspects are taken into account. 

In the following the derivation of these equations is 
described. 

To the Hamiltonian Eq. ([!]) the exciton-exciton inter- 
action H' with a contact potential W is added 

^'=^ E ^<<ak 1+q ak 2 - q . (12) 

ki,k 2 ,q 
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The interaction matrix element W between the excitons 
is in general given by a momentum-dependent expecta- 
tion value of the Coulomb interactions between the various 
point charges in the two excitons containing both the at- 
tractive direct and the repulsive exchange interaction. In 
a low-temperature exciton system the small momentum 
transfer dominates. In this limit the direct interaction can 
be neglected and the exchange integral can be calculated 
analytically . This yields a repulsive contact potential 
with the coupling constant 

26 

W = —™lE r . (13) 

Here do and E r are the Bohr radius and the Rydberg en- 
ergy, respectively. Note however, that the possibility of a 
biexciton formation is lost in this long- wavelength approx- 
imation. 

In the HFB approximation all possible contractions 
have to be taken which conserve the one-particle struc- 
ture of the Hamiltonian and its translational invariance. 
Therefore we have to introduce additionally to p and 
the anomalous function Fk(t) — (afc(t)a_/ £ (t)). From the 
viewpoint of the underlying e-h picture, the bosonic HFB 
approximation includes also the dynamics of four-particle 
correlations (two electron and two hole operators). There- 
fore, correlations beyond the usual electronic HF approx- 
imation are taken into account. 

To derive the kinetics, the Heisenberg equations are 
iterated to first order in W and to second order in g. 
The mean-values of higher orders are then factorized as 
explained in section ||, only that now all possible con- 
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tractions are taken into account including the anomalous 
function F k . The phonon part (collision terms) is treated in 
the Markov approximation and contributions from prin- 
cipal value integrals have been neglected. Due to the fact 
that the number of possible second-order phonon contri- 
butions is quite large, all anomalous functions F k are ne- 
glected in the collision terms. We will show later that the 
contributions of the anomalous functions to the collision 
terms within the same approximation scheme make the 
condensate unstable. In the language of diagram theory, 
our approximation includes the HFB diagrams |l4| and a 
selection of second-order phonon scattering diagrams. 
The resulting equations are 



d_ 

di 



W 
1 i 

o . f^k \coll 

at 



V 



(14) 



-F k = -1 

dt n 



2e k + 2W{n + ^J2 n i^ 



1)W 



(15) 



dt P h 



w . 

eo + -7r{n 



V 



lk Wp *v^ F * + lk dE+ m p ^ 



(16) 



^n k \ coU = ^2al-q{S(e q -e k - 7kJ k _ q ) 

q 

x [N k _ q n k (l + n q ) - [N k - q + 1)(1 + n k )n q ] - (k ^ q)} 
+ —glS(e k - huj k )no{N k - n k ) (17) 



-Q t p\coll = -—p^g 2 q 8{e q -Tiuj q ){N q ~n q ) (18) 
q 

g- t F k \ coU = ~-^g 2 k 5{e k - fiu k ) [p 2 (N k ~ n k )] (19) 



First we demonstrate that the addition of scattering 
terms proportional to F k in the collision terms of F k Eq. ( |l9| ) 
destroy the condensate, while omitting them leads to a sta- 
ble condensate solution. In Fig. [5] a solution of the equa- 
tions is shown above the critical density. While the coher- 
ent density no = \p\ 2 goes to zero if the collision terms pro- 
portional to F k (F 1 -terms) are included (see curves 3,4), 
it condenses if only terms without F k (F -terms) are con- 
sidered as done in Eq.(|19|) (see curves 1,2). For subcritical 
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Fig. 5. Density n and no versus time t above the critical den- 
sity with n tota i — 8.3 ■ 10 16 cm -3 and T p h — 1.5K. Curve 1 
and 2: F neglected(F°) in collision terms. Curve 3 and 4: first 
order contributions(F ) included in the collision terms. Hori- 
zontal lines are n c = 5.2 • 10 16 and ntotai —n c — 3.1 • 10 16 cm' 3 
respectively. 

excitation the condensate decays in both approximations. 

The reason for the numerically obtained condensate 
destroying properties of the anomalous contributions in 
the collision terms are not yet fully understood. It can 
be easily seen analytically that these terms are not in ac- 
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cordance with the condensate solutions of the free Bose- 
gas. However they are obviously also not yet fully com- 
patible with the condensation-solution of the interacting 
Bose-gas. More refined approximations like the inclusion 
of energy renormalizations in the collision integrals are 
needed. An argument for omitting the contributions of the 
anomalous functions to the second order scattering terms 
has been proposed in Ref. [^o). Here it is argued that an 
adiabatic elimination of Fk yields for a scattering term 
containing one anomalous function a term of third-order 
oc Wg^ , which is inconsistent with the considered second- 
order scattering. Although this argument is not rigorous 
in a self-consistent theory, all collision terms proportional 
to Fk are disregarded in the following. 

We compare the prediction for the kinetics of the pure 
phonon model Eqs.@,(§) and the full model Eqs.@jl9|) 
in two figures once below (Fig. |6|) and once above (Fig. ^) 
the critical density. 




with HF 



without HF 



n>n c 



-50 50 100 150 200 250 300 350 400 
time(ns) 

Fig. 7. Excitation at n > n c , comparison HF, without HF 

The results of the two models concerning dephasing 
and relaxation can hardly be distinguished. We conclude 
that all our results in the pure phonon model are stable 
against Hartree-Fock-Bogoliubov correlations. 
Additionally the introduced anomalous function Fk shows 
a threshold behaviour as expected in a phase transition. 
This is shown in Fig. || where we plotted J2 k \Fk(t = oo)| 
versus the excitation density. 
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density (10 cm ) 



50 100 150 200 250 300 350 400 
time ( ns ) 

Fig. 6. Excitation at n < n c , comparison HF, without HF 



Fig. 8. Asymptotic sum of the anomalous function Fk 

To conclude this chapter, we presented a model which 
does show the well-understood dephasing and condensa- 



O.M. Schmitt et al.: Dephasing of excitons in the vicinity of n c 



tion properties of the earlier introduced simpler boson 
model, but additionally can be applied to experiments 
where the coherence plays a more important role than 
in simple studies of the dephasing. 

7 Conclusion 

It was demonstrated that the dephasing of an exciton po- 
larization due to phonon scattering can have surprising 
features, when the BEC is approached. For subcritical ex- 
citation with a laser pulse, the polarization dephasing rate 
slows down with increasing density in the vicinity of the 
BEC with a quadratic power law. For supercritical exci- 
tation the expected exponential dephasing changes to a 
power law relaxation approaching a finite value. We hope 
that these results may stimulate corresponding experi- 
mental investigations. In the last chapter it was shown 
that condensation kinetics is very sensitive to the choice 
of specific scattering terms (diagrams). A model is con- 
structed which includes both the incoherent scattering due 
to phonons and the HFB correlations. It has the same in- 
coherent properties as the pure phonon model, but may be 
also applied to experiments, where coherent correlations 
are important. 

This work has been supported by the DFG in the framework of 
the Schwerpunktprogramm Quantenkohcirenz in Halbleitern. 
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